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Abstract 

One limitation on the performance of optical traps is the noise inherently present in every setup. 
Therefore, it is the desire of most experimentalists to minimize and possibly eliminate noise from 
their optical trapping experiments. A step in this direction is to quantify the actual noise in the 
system and to evaluate how much each particular component contributes to the overall noise. For 
this purpose we present Allan variance analysis as a straightforward method. In particular, it allows 
for judging the impact of drift which gives rise to low- frequency noise, which is extremely difficult to 
pinpoint by other methods. We show how to determine the optimal sampling time for calibration, 
the optimal number of data points for a desired experiment, and we provide measurements of how 
much accuracy is gained by acquiring additional data points. Allan variances of both micrometer- 
sized spheres and asymmetric nanometer-sized rods are considered. 
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INTRODUCTION 



Optical tweezers are the perfect nano-tool for single-molecule manipulation and investi- 
gations [1]. With a correctly chosen wavelength they are nearly non-invasive and can be used 
to manipulate entire living microorganisms or track organelles inside a cell [4]. For 

single-molecule investigations, often a handle, for example in the form of a dielectric micron- 
sized object, is attached to the molecule of interest. Optical tweezers can then follow the 
motion of the handle and calculate the forces exerted on the handle with sub-piconewton 
resolution. The handle is often a polystyrene sphere with a diameter of a couple of microm- 
eters, but even nanoparticles such as gold spheres^jsl, gold nanorods [6], spherical silver 
nanoparticles W , or even individual quantum dots [8] can be individually trapped and used 
as force transducers. 

For single-molecule experiments it is extremely important to measure the distances moved 
and forces exerted by the single molecule as accurately as possible [9]. Therefore, it is crucial 
to minimize or eliminate noise and drift and to perform force calibration as accurately as 
possible. Much effort has been put into minimizing noise, for instance entire setups have 
been covered to eliminate pressure fluctuations; the equipment is most often placed on an 
optical table and sometimes even on a foundation which is separated from the rest of the 
building. A successful way to eliminate drift is by using a laser beam parallel to the trapping 
laser to track the motion of a feducial marker which is attached in proximity to the handle 
of interest and thus subject to a similar drift [l^. This setup, however, requires the use of 
at least two laser beams and two independent detection systems. Another way to reduce 
drift that is often employed in optical trapping setups is to move the system of interest away 
from any surface subject to signiflcant drift. This could be done by using a dual-trap setup 
and suspending the molecule of interest between two individual handles. But how eflicient 
are these methods? And which types of noise should one really worry about in optical 
trapping experiments? To answer these questions it is essential to be able to quantify the 
noise introduced by each part of the equipment or surroundings. 

Fourier analysis is an excellent tool to calibrate optical tweezers on-the-fly. Furthermore, 
peaks occuring in the high frequency part of the power spectrum can often easily be traced 
backwards to a particular noise contribution. However, Fourier analysis is not optimal 
for identifying and quantifying low- frequency noise, which drift, as inherently present in 



experiments, produces. However, Fourier analysis has been used to judge the noise stemming 
from particular experimental settings [ll, 12], but for such an analysis various assumptions 



have to be made about the bandwidth of the integration. Also, the regular positional variance 
is often used as a measure for noise. However, the normal variance does not converge for 
purely stochastic types of noise, such as white noise. 

In this Proceeding, we propose Allan variance analysis as a simple and efficient tool to 
pinpoint and quantify noise in optical trapping facilities. The Allan variance of the positions 
visited by an optically trapped particle can be calculated on-the-fly during an experiment and 
used, e.g., to determine the optimal length of a time series for accurate calibration Q- 
By comparison to simulated data with no drift present it has been explicitly shown that 
Allan variance is an excellent tool to pinpoint and quantify low-frequency drift Q,^. Allan 



variance analysis has been specifically used to verify that a CMOS camera had sufficient time 
resolution to track an optically trapped particle [131]. In addition, Allan variance analysis is 
able to reveal noise contributions from commonly used photodiode-based detection systems 
in optical trapping experiments jl^, thus providing a platform for a qualified choice between, 
for example, a position-sensitive rather than a conventional quadrant photodiode for certain 
sets of experiments. Moreover, Allan variance has revealed the impact of the piezo stage, 
the acoustic noise in the laboratory, and the geometry and stability of the sample chamber 
on the noise spectrum. Complementary and additional to Fourier analysis, Allan variance 
analysis provides a basis for an optimal setup design [l^. In this Proceeding, we review 
some of these results. In addition, we provide information about how the Allan variance of 
a measurement can be improved by acquiring additional data points and about the Allan 
variance of an asymmetric gold nanorod. 



METHODS 



By focusing a laser beam into a measurement chamber that is mounted onto a piezo stage 
inside an inverted microscope (Leica DMIRBE), the optical trap is created. For trapping 
polystyrene speres (Bangs Laboratories, diameter (800 ± 10) nm) we used a water immersion 
objective (Leica, HCX, 63x, NA=L20) at a color correction at its lowest setting (0.13 mm) 
and a measurement chamber of two cover slips (bottom: thickness 0.13^0.16 mm, top: 1 mm) 
sandwiched together by double sticky tape. The actual trap was formed in the middle of 




beam 



FIG. 1: Schematic drawing of experiments. A highly focused infrared laser beam traps a 
polystyrene sphere inside a custom-made measurement chamber that is mounted on a piezo stage. 
The forward scattered light is collected onto a photodiode that transmits the sphere's position into 
an acquisition system. 

the chamber (total height (95 ± 5) fim). For trapping gold nanorods an oil immersion 
objective (Leica, HCX PL Apo, 63?, NA =1.32) was used and the immersion oil was chosen 
to compensate spherical abberations flG^. Here, the measurement chamber was custom- 
made from two cover slips (bottom: thickness 0.16-0.19 mm, top: 1 mm) held together by 
parafilm. The nanorods were trapped 5 /xm above the bottom. Their dimensions were for 
the longer axis x = (63.8 ± 7.4) nm and for the shorter axis y = (37.3 ± 5.0) nm (verified 
by transmission electron microscopy prior to measurements). In all cases, measurement 
chambers were sealed with vacuum grease to prevent evaporation. To ensure the optimal 
pointing stability of the laser beam, we switched on the laser at least one hour prior to 
experiments. All measurements were done at room temperature. 

A schematic drawing of the essential parts of the experiment is shown in Figure [H The 
forward scattered light is focused onto a position-sensitive photodiode (Pacific Silicon Sensor, 
DL100-7PCBA3) or onto a quadrant photodiode (Hamamatsu, S5981). Its output in voltage 
is connected through an amplifier, a low-pass filter of 100 kHz, and an acquisition card 
(National Instruments, PCI-6251) to a computer. By utilizing our data-streaming software 



that was custom-made in Labview^ 



17|, we were not limited by the amount of acquired 



data (tested 2 h at 1 kHz, and 22 min at 100 kHz). 

Typically, for trapped spheres, the analysis was applied to adjacent time series of 2^^ 
positions acquired at various acquisition frequencies. For nanorods, the time series consisted 
of 2^^ positions acquired at 22 kHz. All additional filtering was carried out posterior to 
acquisition. For visual observation of the spheres right before and after experiments, a CCD 
camera (Sony, XC-ES50, 25 Hz) was used, whereas the nanorods could not be visualized 
optically. 

Calibration. An optically trapped object experiences a harmonic force F = —nx with 
the trap stiffness k and the distance x from the equilibrium position. Thus, k characterizes 
the thermal motion of the trapped object. For the analysis, the Langevin equation is typi- 
cally solved and Fourier transformed. The result is a positional power spectrum that allows 
for finding the ratio between n and the friction coefficient 7, i.e. the corner frequency fc'. 

I' - h 

In case a sphere is trapped far away from any surface, the Stokes law gives: 7 = 67rrr7, where 
r is the radius of the sphere and t] the viscosity of the surrounding medium; here, in all ex- 
periments, deionized water, r\ = 8.9- 10""^. The nanorods have the shape of sphero-cylinders. 
We approximated their shape as cylinders in order to calculate their drag coefficients ac- 
cordingly. [6] The conversion factor /?, which relates the distance measured in volts to the 
distance travelled in meters by the trapped object, was found by comparing theoretical and 
experimental diffusion constants fisl . 

In a first step, we calculated the Allan variance to obtain the optimal measurement 
time for calibration. Secondly, we calibrated conditionally independent intervals of the 
time series with that particular length. We used the power-spectrum method as described 
previously [19] with the freely available program \2Q\. Finally, we calculated various types 
of variances. 

For typical noise phenomena found in nature, classical variances do not converge. For 
example, purely stochastic noise does not converge for the normal variance: 

a'{T) = l{{x,-x)\, (2) 

where x denotes the mean of the time series and Xi the mean of positions within the interval 
of length r. 



The Allan variance is designed to converge for most naturally occuring noise. [2j| Given 
a time series consisting of elements and a total measurement time of tacq = /acq^, the 
Allan variance is defined as: 

W = \ {{x^+l-Xif)^, (3) 

where Xi is the mean of the measurement interval r. In words, the Allan variance is half of 
the mean of the squared differences of neighboring intervals of a given length. 

One can trade the Allan variance's conditional independence of neighboring intervals to 
gain a much smaller statistical error. For the overlapping Allan variance, one simply 
calculates all possible diflFerences of neighboring intervals in a given time series. For a more 
comprehensive discussion we refer elsewhere jl4|. 

As thermal limit for an object in an spatially confined trap, the standard error of an 
object's position averaged over the time interval r is [1^: 



with ksT being the thermal energy. This limit cannot be bettered by any measurement of 
an object trapped by a single beam. Nevertheless, for so-called dual-beam traps this limit 



is about a factor 1.19 smaller 



13|. 



Allan variances were calculated with a custom-made Matlab program [22[]. We measured 
Allan variances for various objects and parameters (trap stiflFness, acquisition frequency, 
piezo) by acquiring three time series. Then the parameter was altered and we repeated the 
measurements. If not stated otherwise, we plotted the overlapping Allan variance in a log- 
log plot. All results stated in the following section are explicit and reproducible, although 
we chose to plot only particular data sets in order to keep a clear representation. 



RESULTS AND DISCUSSION 



In the majority of the experiments a polystyrene sphere was trapped in an aqueous 
environment in the center of the measurement chamber such that the distance to any surface 
was significantly larger than the diameter of the sphere. Using the position sensitive diode 
the time series of the sphere was recorded and used for a calibration routine which returned 
the corner frequency /c, the trap stiflFness and the conversion factor (3. Figure [2] (a) shows 
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FIG. 2: Time series analysis of an optically trapped polystyrene sphere, k = 67.7 /xm/pN. (a) 
Position as a function of time, (b) Histogram of positions, fitted by a Gaussian distribution, (c) 
Power spectrum of positions. The full line is a fit to the grey region which incorporated aliasing 
and filtering effects. 

the time series of the positions visited by a trapped particle. The histogram of all positions 
visited is plotted in (b), and is very well- fitted by a Gaussian distribution. Figure [2] (c) 
shows the corresponding power spectrum. Here, the light grey part is fitted by a Lorentzian 
function (black full line) which takes into account the filtering effect of the photodiode [23] 
and aliasing using programs described above [2^. Though the system is subject to 
low-frequency drift, this does not show neither in the time series, nor in the histogram, nor 
in the power spectrum, not even when compared to a simulation of the situation without 
drift using similar physical parameters [14]. 

From time series as in Figure [2] the Allan variance can be calculated using Equa- 
tion ([3j) [22(] . A typical result is shown in Figure [3] where the Allan variance is plotted 
as a function of data acquisition time for a strong, /^=63.9 pN//xm (black), and a weak, 
/^=34.7 pN//xm (grey), optical trap. The thicker lines denote the overlapping Allan variance 
for the same time series. The two dashed lines with slopes of —1/2 correspond to the thermal 
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FIG. 3: Variances of a polystyrene sphere trapped strongly, hc—63.9 pN//xm (black), or weakly, 
K=34:.7 pN//xm (grey). The lighter graphs denote the Allan variances, the thicker ones the over- 
lapping Allan variances. Dashed lines are the thermal limits, dotted lines the normal variances. 

limits (Equation H]). For short measurement intervals the Allan variance is smaller than the 
thermal limit, this is due to the correlation of the data points. The maximum of the Allan 
variance is at ttTc where 

.r. = ±-. (5) 

For measurement intervals where the Allan variance is larger than the thermal limit, the 
data points are not correlated. The observation that the Allan variance is very close to 
the thermal limit in the measurement time span between 100 ms and 1 s is evidence of an 
extremely stable setup. The Allan variance has a global minimum at around 1^10 s, the 
exact post ion of this minimum is dependent on the trap stiffness. This global minimum 
denotes the optimal measurement time for calibrations. This minimum is the time where 
the Gaussian distributed parameters have been measured for a time long enough to allow for 
their fairly precise determination while the time interval is still short enough that drift is not 
yet a significant problem. Noticeable also, is the fact that the stronger trap has a lower Allan 
variance, hence, a better accuracy than the weaker trap. The dotted lines in Figure [3] denote 
the normal variance (Equation [2D. The normal variance does not converge and consequently 
does not provide the necessary resolution to determine the optimal measurement time. 
Figure [H (a) shows the Allan variances as a function of the number of acquired data 
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FIG. 4: Frequency dependence of Allan variance and its negative derivative for a sphere trapped 
with pN//im. (a) Allan variance as a function of number of acquired data points. The full 

line denotes the limit of how much positional information a single data point could possibly hold, 
(b) Negative derivative of Allan variance with respect to the number of data points. 

points. A polystyrene sphere is trapped with tv=67.7 pN//xm while the acquisition frequency 
is parameterized. The tested frequencies range from 10 Hz to 100 kHz. For all frequencies, 
the acquisition of additional data points does not increase the accuracy above a certain 
threshold, for instance at 2^^ for 10 kHz, because the Allan variance shows an absolute 
minimum. Furthermore, the closer a graph stays to limit of maximum information per 
individual data point, i.e. quasi the 'thermal limit' in this way of representation, the more 
valuable is the acquisition of an additional data point. Figured] (b) shows the negative values 
of the derivatives of the Allan variances with respect to the number of data points - it gives 
the increase in accuracy by acquiring more data points. For a specific acquisition frequency, 
it passes through zero when its sampling has been optimal and for more points drift would 
start to dominate. Consistently, these graphs show that if the sampling frequency is 10 kHz, 
essentially nothing is gained by acquiring more than 2^^ data points. 
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FIG. 5: Contribution from the piezo stage to the Allan variance. The full line shows the Allan 
variance of a trapped sphere when the piezo stage is off, the dashed line is the Allan variance when 
the piezo is switched on. The dotted line is the Allan variance of the position output from the 
piezo control box with a significant noise contribution over a broad frequency band indicated by 
grey shading. 

Often, the measurement chamber of an optical trapping experiment is mounted on a 
piezo-electric stage. To investigate the possible contribution from the piezo stage to the 
noise spectrum we measured the Allan variance for two similar trapping experiments. In 
Figure [5] one sees the first case with the piezo switched off (full black line) , and the second 
case with the piezo switched on (dashed line). There is a distinct peak in the Allan variance 
at around 10 s which seems to originate from the piezo stage. The lower dotted line in 
Figure [5] is the Allan variance of the position of the piezo stage as given by the piezo control 
box. This noise spectrum also peaks at around 10 s and hence supports the conclusion that 
the piezo stage itself contributes to the noise spectrum over a broader low-frequency band as 
indicated by the grey shading. At short measurement times the Allan variance of the piezo 
output has some oscillations which correspond to the odd-numbered divisors of the piezo's 
resonance frequency of 100 Hz. The fact that the piezo contributes to the noise spectrum 
implies that if it is not strictly needed for a particular experiment, it should be switched off. 

So far, all figures and results presented both in the present Proceeding and in litera- 
ture [3, Q have addressed optical tweezing of micrometer-sized spherical objects. It is. 



however, also possible to trap significantly smaller and non-spherical objects, such as gold 
nanorods. Gold nanorods as thin as 8 nm and with aspect ratios up to 5.6 have been 
optically trapped [6]. These asymmetric nanorods align inside the trap with their longest 
direction along the electrical field vector of the trapping laser and the spring constant corre- 
lates directly with the polarizability of the rod. Figure [6] shows the Allan variance calculated 
from an experiment where an individual gold nanorod (long axis (63.8 ± 7.4) nm, short axis 
(37.3 ± 5.0) nm) was optically trapped. The black trace is the Allan variance along the 
longest dimension of the rod, the grey trace along the shortest dimension of the rod. The 
thermal limit for the shortest dimension is lower as we consider here r to be the short 
axis (Equation (111). Figure [6] also points out that for long measurement times, the Allan 
variance is smaller for a stronger trapping stiffness. Moreover, it shows that the optimal 
measurement interval is shorter for lower thermal limits rather than for the trap stiffness. 
The optimal measurement intervals are on the order of a tenth part of a second, which is 
significantly shorter than the optimal measurement times for the much larger polystyrene 
spheres, typically on the order of seconds (Figure [3]). At measurement times longer than 
the absolute minimum, there are distinct peaks in the Allan variances of both directions. 
As the measurement chamber was as narrow as 30 mm x 5 mm, it might be possible that 
those peaks originate from the geometry of the chamber [14]. The dotted lines denote the 
normal variances of the plots, which, also here, do not reveal any information about, e.g., 
optimal measurement time or the impact of drift. 



CONCLUSION 



We presented Allan variance analysis as an excellent tool to quantify noise in experi- 
ments where micro- and nanometer-sized particles were optically trapped. Furthermore, we 
pointed at our proper software solutions that can easily be used for on-the-fiy analyses to de- 
termine important measurement parameters, as for instance the optimal measurement time 
for calibrations, or the optimal number of acquired data points for a particular experimental 
setting. 

Fourier analysis can conveniently be used to perform calibration procedures. Further- 
more, its use to pinpoint high-frequency noise is excellent. However, for low-frequency noise 
which is inherently present in every experiment, Allan variance analysis is superior. This 




FIG. 6: Variances of optically trapped gold nanorods. The black curve is along the longest direction 
of the rod, k = 2.3 pN//im, the grey curve along the shortest direction of the rod, hi = 2.1 pN//im. 
The dotted lines are the corresponding normal variances, the dashed lines the thermal limits. 

method has been used to quantify the noise contribution from photodiode-based detection 
systems, the influence of the noise on chamber stabihty and geometry 1^, and, as reviewed 
here, the noise contribution from the piezo stage. Allan variance analysis was mostly used on 
positional time series from spherical, micron-sized polystyrene spheres, but it can be utilized 
even for nanometer-sized gold rods. Here, the overall noise contribution was asymmetric, 
strongly correlated to the alignment of the nanorods inside the trap. 

We think that Allan variance analysis, complementary and additional to Fourier analysis, 
allows for the development of a common standard in research that accesses optical tweezers, 
as it contains the possibility to compare noise and drift in different experiments, settings, 
setups, and even laboratories. 
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